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We consider heat conduction in a ID dynamical channel. The channel consists of a group of 
nonintcracting particles, which move between two heat baths according to some dynamical process. 
We show that the essential thermodynamic properties of the heat channel can be evaluated from the 
diffusion properties of the underlying particles. Emphasis is put on the conduction under anomalous 
diffusion conditions. 
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The link between thermodynamic phenomena and mi- 
croscopic dynamical chaos has been a subject of interest 
for a long time [1]. Examples of such a relationship is 
deterministic diffusion [2], where local dynamical prop- 
erties, such as stability (or instability) of several fixed 
points can change the global diffusion transport from 
normal to anomalous [3]. Another question which has 
attracted a lot of attention recently, is the problem of 
heat conductivity in deterministic extended dynamical 
systems [4-8]. A large number of models have been pro- 
posed in order to understand the conditions under which 
a system obeys the Fourier heat conduction law [4] . Re- 
cently, a new class of ID models, "billiard gas channels", 
has been proposed [5-8] . These channels consist of two 
parallel walls with a series of scatterers, distributed along 
walls, and noninteracting particles that move inside. The 
two ends of the channel are in contact with the heat 
baths. By changing the shapes and positions of scatterers 
it is possible to change the conductivity of the channel 
[5-8]. 

In this Letter we show that such billiard gas channels 
belong to a wider class of models, which we call dynam- 
ical heat channels. The absence of interactions between 
particles and the independence of the particle dynamics 
on kinetic energy allow a complete separation between 
the thermodynamic aspect, which is governed by proper- 
tics of the thermostats, and the dynamics inside channel, 
which is governed by diffusion properties within the chan- 
nel. All the essential information on heat conductivity of 
such a dynamical heat conductor can be obtained from 
the diffusion properties of channel. 

Dynamical heat channel. To model dynamics within 
the channel we consider N particles that move along di- 
rection X, following the dynamical equations of motion 

x = f(x,i), A ex, (1) 

where the function f can be either deterministic or ran- 
dom. 

To consider transport of heat, two heat baths with tem- 
peratures T + and T_ are attached to the left and right 
ends of the channel. The heat bath is characterized by 
a velocity probability density function (pdf), PrivTh), 
where vtk is the thermal velocity. After colliding with 



the heat bath, the particle is ejected back to the channel 
with a velocity vtk which is chosen from Pt{vtk)- 

As particles do not interact, the dynamics of the en- 
semble inside the channel can be described by a long 
trajectory of a single particle and the flux should be 
rescaled by the factor N. The trajectory of a particle 
is independent of the particle velocity vtk- The only dif- 
ference between "hot" and "cold" particles is that the 
hot ones cover the same trajectories faster then the cold 
ones. The velocity vj-h does not change during the prop- 
agation through the channel and it can be interpreted as 
a temperature "label" of a particle. The dependence of 
the dynamics in Eq. (1) on vtk can be taken into account 
by introducing a scaling factor for the time t — > t/vxh- 

Due to the separation of the thermodynamic charac- 
teristics from the dynamical ones, the proposed approach 
is not limited to Hamiltonian systems only [5-8] . We only 
assume that the dynamics inside the channel has a diffu- 
sional character, and can be characterized by evolution 
of the mean square displacement (msd) 

(X 2 (t))~tf\ (2) 

This diffusion can be normal (a = 1), subdiffusive (a < 
1) or superdiffusive (a > 1) [9]. 

Following Ref. [5], heat transfer by a particle through 
the channel is 

M(t) M(t) 

Q(t) = J2 A ^ = £ *i • i.Ef - Ep% (3) 
j=i j=i 

where E™ and E° ut are the energies before and after the 
j-collision with the heat bath, qj is the direction factor, 
qj = 1 if the (j — l)-collision is with the hot end, and 
qj = —I in the case of the cold end. M(t) is the total 
number of collision events during time t. In the case of 
normal heat conductivity, Q grows linearly with t and 
the heat flux is defined by 

J = Um t ^ 00 9¥L. (4) 

Let us start from a situation where the particle is ini- 
tially located at the hot end. During diffusion it can 
come back and collide with the hot bath again. Accord- 
ing to Eq.(3) this event, on average, does not lead to 
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heat transfer. But when the particle reaches the opposite 
cold end Q increases on average as / °° ^^[Pr+ivTh) — 
P T (vTh)]dvTh- After that the process is reiterated start- 
ing from the cold end. Thus, the problem of heat transfer 
is reduced to the problem of diffusion in a finite inter- 
val under reflecting and absorbing boundary conditions. 
As the initial condition we assume that at t = the 
particle is located at the reflecting end. The average 
time r needed to reach an absorbing boundary is the 
first moment of the pdf <j)(f) for first arrival times. To 
take into account the effect of thermodynamic velocity 
on the time r should be rescaled, as mentioned above, 
t — > t J °° -^-PT ± (vTh)dvTh ■ Due to the absence of 
mass flux, the number of transitions from left to right 
and vise versa should be the same. Finally, we obtain the 
following equation for the one-particle heat flux through 
the channel of length L 



J(L) 



1 Jo v Th \ P T+ {vTh) - PT- (vTh)}dv T h 
/o°° vth [ Pt + ( Wr '») + Pt - ( V Th)]dv T h 



■ (5) 



For an ensemble of particles the heat flux should be writ- 
ten as J e ns — N ■ J(L). Here we take into account that in 
order to keep a density of particles in the channel fixed, 
-/V should be proportional to L; namely N cx L. 

Equation (5) demonstrates a complete separation be- 
tween thermodynamic and dynamical aspects. Without 
loss of generality we consider the "delta" -heat bath with 
a simple pdf, P T (v Th ) = S(v Th - \/2T) [8]. 

As a model for the dynamics in the channel we con- 
sider continuous time random walks (CTRW) [10] in a 
discrete lattice. This allows to cover the spectrum of dif- 
fusion regimes, from subdiffusion to superdiffusion, and 
to describe kinetics of deterministic Hamiltonian [11,12] 
and dissipative [3] systems. 

The CTRW-model is a stochastic process which rep- 
resents an alternating sequence of waiting and jumping 
events [10]. A particle waits at each point for a time cho- 
sen from a waiting time pdf tp w (t), and makes a jump 
(flight) to the left or right with equal probabilities. The 
jumps are characterized by the pdf ipf(x,t), the proba- 
bility density to move a distance x at time tin a a single 
fligth event. We consider power laws for both pdf's[10] 



ip w (i) <~ t~ 



i> f (x,t) 



~ x 



-ii- 



(6) 



Depending on asymptotic properties of ip w (t) and 
i>f(x, t), we have regimes of superdiffusion, normal diffu- 
sion or subdiffusion [10]. 

Superdiffusion. To achieve superdiffusion with a finite 
msd (x 2 (t)) there should be a correlation between the 
length and duration of the individual flights. Such corre- 
lation leads to the model of the Levy walks with spatio- 
temporal pdf [10] 
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FIG. 1: pdf 4>(t) for first arrival times for (a) superdiffusion 
with 7 m = 2.6, 7/ = 1.6 for L — 100. Inset showns the 
scaling of the exponent s with the channel length L, L~ lf . (b) 
Subdiffusion, j w = 0.6, 7/ = 2.6. Straight lines correspond 
to the asymptotic behaviour (see text). 



that corresponds to flights with a constant velocity v. 
Here we assume that j w > 1, so there is a finite mean 
waiting time. Depending on 7/ we distinguish among 
three regimes of diffusion [10], according to the exponent 
a in Eq.(2), 



< 7/ < 1 
" = : < ■'> - If, 1 < 7/ < 2 , 
2<7/ 



(8) 



V/Or, t) = iP f (t)6(\x\ - vt), ^ f (t) - t-ir 



(7) 



In order to derive the dependence of the heat flux J 
on the length L of the channel we start from a consider- 
ation of the survival probability <&{t) for a particle walk- 
ing on the finite interval of the length L bounded by 
reflecting and absorbing boundaries. We find that for 
Levy walks the survival probability is cx e~ st , where 
s cx L^ 1 for the ballistic motion (a — 2), s cx L~ 7/ 
for superdiffusion (1 < a < 2) , and s cx L~ 2 for nor- 
mal diffusion (a = l)(sce inset in Figla). Then, the 
pdf for first arrival times, <j){t), is <j>(t) <~ — &(t) cx e~ st 
(Fig. la) and has a finite first moment, r. As a result, in 
asymptotic limit the heat, Q, grows linearly with time 
t, Q{t) - (AE/r)t (line (3) in Fig.2), where the mean 
arrival time is r = J°° t<p(t)dt cx s _1 (Fig. 3a). Thus, in 
the case of superdiffusion, due to Levy walks, we get a 
normal linear in time heat conductivity. 

Finally, for the one-particle heat flux, following Eqs. 
(5) and (8), we arrive at the following equation for J(L) 
in term of the msd exponent a, 

( 1, a = 2 

J(L)<xL- t3 , p=l 1} = Z- a , l<oc<2, (9) 
[2, a=l 

In Fig3.b we show the numerical results for J(L) obtained 
for 7/ = 1.6 (stars) and 7/ = 2.4 (squares). The results 
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FIG. 2: Evolution Q(t) vs t for subdiffusion, 7^ = 0.6, 7/ = 
2.6 (curve (1)), normal diffusion, 7™ = 1.6, 7/ = 2.4 (curve 
(2)), and superdiffusion, 7™ = 1.6, 7/ = 1.6 (curve (3)). 
Curve (4) corresponds to rescaled, for one particle, Q(t) for 
an ensemble N = 100 in the case of competition between 
anomalous flights and localized events, 7™ = 0.4, 7/ = 1.2. 
Straight lines correspond to asymptotics (see text) For all 
cases the length of channel is L = 1000, T+ = 2 and T_ = 1. 



are in good agreement with the scaling law suggested by 
Eq.(9) (straight lines in Fig. 3b). 

Then taking Eq.(9) into consideration the thermal con- 
ductivity k — —J ens (N)/VT, shows the following behav- 
ior for 1 < a < 2, 

kKl^'al"- 1 . (10) 

This diverges as one goes to the thermodynamic limit 
L — > 00. For the ballistic case a — 2, k oc L and for 
a = 1 k oc const as expected. 

The Levy walk model is an adequate approach for mod- 
cling of Hamiltonian kinetic in a mixed phase space [12]. 
Previous numerical results which have been obtained for 
Hamiltonian billiard channels, show that a = 1.3 corre- 
sponds to the flux exponent (i = 1.72 [7], and a = 1.8 
to (3 = 1.178 [8], which are in good agreement with the 
relation in Eq.(9). 

In the case of Levy flights [10], unlike the Levy walks, 
the velocity is not introduced explicitly. Here we assume 
that all flights have the same duration tf and use the 
following decoupled representation for flights pdf, 

1>f{x,t)=1> f (x)6{t-tf), 4> f (x)~x-<'-\ (11) 

In this case the msd diverges and one is tempted to use 
(|a;(i)|) 2 to characterize the dynamics. In the asymptotic 
regime this quantity scales as (\x(t)\) 2 <~ i 2 / 7 ^[10]. As in 
the previous case, here the exponent 7/ determines the 
dynamics in the channel and the dependence of flux J(L) 
on L. The scaling of the mean arrival time, r cx L ls , is 
consistent with the result obtained from a direct solution 
of a fractional Fokkcr-Planck equation [14] . We note that 
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FIG. 3: (a) The mean first arrival time r versus L; (b) one- 
particle flux J(L), Eq. (5), versus L ( T+ = 2 and T_ = 1). 
Results correspond to different diffusion regimes with subd- 
iffusion, 7m =0.6, 7/ =2.6 (circles), normal diffusion j w = 
1.6, 7/ = 2.4 (squares), superdiffusion j w = 1.6, 7/ = 1.6 
(stars), and Levy flights, 7™ = 1.6, 7/ = 1.6 (triangles). 
Straight lines correspond to asymptotics (see text). 

using of (|ir(i)|) 2 instead of msd gives a different, yet un- 

— -3 

physical scaling, J(L) a L 1 / , with an exponent which 
is close to that in the scaling in Eq. (9) and coincides 
with it at the points 7/ = 2 (normal diffusion) and 7/ = 2 
(ballistic diffusion). 

Subdiffusion. In the case of subdiffusion 7^ < 1 and 
7/ > 2 in Eq.(6) and the mean waiting time t w diverges 
[10]. The motion is characterized by long localized events 
and the msd grows sublinearly, (x 2 (t)} ~ i 7 -. The pdf of 
the first arrival time, <fi(t), has a power-law asymptotics 
<f>(t) oc i^ 7 ™" 1 (Fig. lb). This finding is in line with a re- 
sult for subdiffusion under a bias in a seminfinite interval 
with an absorbing boundary [15]. The anomalous char- 
acter of the pdf 4>(t) leads to a divergence of r. Thus, 
in the case of subdiffusion the heat Q{t) carried through 
the channel grows sublinearly with time (see line (1) in 
Fig.2), 

Q{t)~t^. (12) 

Within a traditional definition of the flux, Eq.(4), the 
subdiffusive anomalous heat conductivity can not be dis- 
tinguished from a heat isolator. We note that subdiffu- 
sion regimes can not be achieved in the case of channels 
with Hamiltonian dynamics due to a finiteness of recur- 
rence time [9]. 

Coexistence of long localized events and fligths. In 
the case of competition between fligths and localization 
events the msd exponent a is given by a relation [10] 

a = 2 + lw - lf , < lw < 1, 1< 7/ < 2. (13) 

The pdf for the arrival time, <f>(t), is governed, like in 
the case of simple subdiffusion, by the pdf for waiting 
times, ip w (t). The divergence of mean arrival time r can 
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FIG. 4: Evolution Q(t) vs t for the map in Eq. (14) with 
7™ = 0.5, 7/ = 1.5 (L = 1000, dt = 10~ 4 , T+ = 2 and T_ = 
1). Straight lines correspond to power law dependence in 
Eq. (12). Inset shown numerically obtained relation between 
scaling exponent (3 for one-particle flux and exponent a (7™ = 
1.6 for all cases). Straight lines correspond to f3 — 3 — a. 



The values of the flux exponents /? obtained for the 
map in Eq. (14) are in agreement with the prediction 
in Eq.(9) for superdiffusion (see inset in Fig. 4). In Fig4. 
we show the time evolution of Q(t) for the case of sub- 
diffusion, "f w = 0.5, 7/ = 1.5. The steps in the Q(t) de- 
pendence correspond to anomalously long waiting events 
when particle is trapped near marginal stable fixed points 
x = j + ±,j = 0,....,L [3]. 

In summary we have introduced a class of dynami- 
cal heat conductors. This class includes as particular 
cases recently proposed Hamiltonian billiard channels [5- 
8]. In the absence of interactions between the particles 
and the independence of the dynamics on particle energy 
the proposed aproach goes beyond Hamiltonian dynam- 
ics and allows to express heat conductivity in terms of 
channel diffusion properties. The Hamiltonian character 
of dynamics becomes essential when interactions between 
particles are introduced (like in the case of dynamical lat- 
tices [4]) or a dependence of dynamics on particle energy 
is included. 

Financial support from the Israel Science Foundation, 
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lead to anomalous heat conductivity Eq. (12), even in 
the case of superdiffusion spreading, (x 2 (t)) <~ t a , a > 1, 
Eq. (13). Fig.l (line (4)) shows the dependence Q(t), for 
the ensemble of N = 100 particles (7^ = 0.4,7/ = 1.2, 
a = 1.2), rescaled for one particle. 

One can go continuously among the various diffusional 
behaviors by tunning parameters in an iterated map. As 
an example of deterministic channel we consider a com- 
bined map [3]. This one-dimensional map generates in- 
termittent chaotic motion with coexisting localized and 
ballistic motion events. The map is defined for one unit 
cell by the iterative rules, 

X n +i = ,f(X n ), 
( X + aX'-l, 0<A<i, 

■ f{x) -\x-~a(±-xy, \<x<\. {u) 

The parameters are a — 4 z z/(z + z) and a — A z z/{z + 
z). The localized and ballistic phases are characterized 
by the pdf's in Eq. (6), with the exponents 7/ = (z— 1) _1 
and j w = (z — 1) _1 . Variation of the map parameters z 
and z allows to cover the diffusional regime in the range 
< a < 2 . 

The thermodynamic velocity, vm , can be included into 
the map description. We fix a time step dt and after one 
iteration of the map, Eq.(14), stretch the system time 
following t = t + dt/vTh- The length of channel is deter- 
mined by a number of unit cells, > [X] > L. When the 
particle reaches the boundary cells, [X] = or [X] = L, 
it is randomly placed into the interval 0<A<l/2or 
L — 1/2 < X < L) respectively. This corresponds to a 
diffusive reflection of particle back to the channel. 
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